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Abstract 

We report the first results of relativistic correlation calculation of the spectroscopic properties 
for the ground state of E112H and its cation in which spin-orbit interaction is taken into account 
non-perturbatively. Studying the properties of E112 (eka-Hg) is required for chemical identification 
of its long-lived isotope, 283 112. It is shown that appropriate accounting for spin-orbit effects leads 
to dramatic impact on the properties of E112H whereas they are not so important for E112H + . The 
calculated equilibrium distance, #^=1.662 A, in E112H is notably smaller than i?f pt =(1.738 ± 
0.003) A and i?f lc =1.738A in HgH, whereas the dissociation energy, Z^ alc =0.42 eV, in E112H is 
close to Df pt =0.46 eV and L>= alc =0.41 eV in HgH. These data are quite different from i?^ H =1.829 A 
and -D^ H =0.06 eV obtained for E112H within the scalar-relativistic Douglas-Kroll approximation 
[Nakajima and Hirao, Chem. Phys. Lett., 329, 511 (2000)]. Our results indicate that E112 should 
not be expected to be "more inert" than Hg in opposite to the results by other authors. 
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Introduction. The superheavy element 112 (eka-Hg) was discovered at GSI (Darmstadt) 
in 1996 within the "cold" fusion reaction 1 . The recent observation at FLNR (Dubna) of the 
"a-SF" chain, attributed to 4 sec a-decay branch of 283 112 followed by a 0.2 sec spontaneous 
fission of 279 110 (Ds)2, brought up the question, what species was observed in the previous 
"hot" fusion FLNR experiment 3 -. Moreover, the production of 283 112 in the reaction of 48 Ca 
and 238 U was not confirmed at LBNL (Berkeley) 4 ^. However, a very specific decay mode 
of the short "a-SF" chain offered a unique chance to unambiguously identify 283 112 in a 
chemical experiment. 

To our knowledge, starting from the papers of Pitzei^ and Fricke^ in 1975 it was main! 
suggested by other authors that El 12 behaves rather like a rare gas than Hg. In Ref. 
the confusing conclusions about both relative inertness of El 12 as compared to Hg and 
similarity of E112 with Hg were made in the abstract and conclusion, respectively. The first 
attempt to identify El 12 chemically was made at FLNR^* 1 ^ but no spontaneous fissions were 
detected. It was interpreted as indication of the Rn-like behavior of El 12 as well. 

The chemical experiments on studying properties of El 12 are currently under way at 
FLNR—, similar work is in progress at GSI and PSI (Villigen) 1 ^ that involves the attempt 
to clarify the recent observation of the decay chains and fission products associated with the 
production of E114 and E116 2 -: being their decay product, E112 should be detectable in gas- 



phase chromatograp 
properties (see Refs. 



lie experiments. The experimental study of superheavy element (SHE) 



l3fl3l and references) is very difficult because of their short half-lives and 



extremely small quantities, only single atoms are available for research. In this connection, 
reliable theoretical prediction of their properties based on benchmark ab initio calculations 
is highly desirable. As a first step in the extensive study of chemistry of E112, bonding 
in simple diatomic molecules such as E112H should be studied. The earlier studies of eka- 
mercury fluorides^ 4 " are not so sensitive because both mercury and xenon are known to burn 
in fluorine atmosphere. On the other hand, the Hg2, Xe2 and El 122 dimers^* 1 ^ are Van 
der Waals systems with a small dissociation energy. By contrast, the ground state RnH 
and XeH molecules are not observed in the gas phase, whereas HgH can be obtained by 
radiofrequency discharge in hydrogen and metal vapor (see, e.g., Ref. Ilo). 

It was shown in Ref.[l7|on example of E112 and other SHE's that the errors in calculations 
due to employing the point nucleus (instead of the realistic Fermi nuclear model) reach 0.4 eV 
for transition energies between low-lying states, whereas neglecting the Breit effects leads to 



2 



the errors up to 0.1 eV. The generalized relativistic effective core potentials (GRECP's)^ 
were generated for El 12 and other SHE's 17 which allow one to simulate Breit interaction 
and Fermi nuclear model by economic way but with very high precision^. The accuracy of 
these GRECP's and of the RECP's of other groups was estimated in atomic finite-difference 
SCF calculations with Coulomb two-electron interaction and point nucleus as compared to 
the corresponding all-electron Dirac-Fock-Breit calculations with the Fermi nuclear model. 
It was justified and checked in Refs. 17 3 that the GRECP method allows one to carry out 
reliable calculations of SHE's and their compounds within the level of "chemical accuracy" 
(1 kcal/mol, 0.043 eV, or 350 cm -1 for valence transition energies) when the valence and 
outer core shells are appropriately treated. Hence, the overall accuracy of calculations 
in heavy-atom molecules is limited, in practice, by current possibilities of the correlation 
methods and codes and not by the GRECP approximation. 

In this paper, we present GRECP calculations of spectroscopic constants for the ground 
states of the E112H molecule and its cation exploiting molgep code^. To our knowl- 
edge, only three calculations on E112H and E112H + were published ^ 14121 . In Ref. bjj, the 
third order Douglas-Kroll (DK3) method was applied to calculation of E112H and its ions. 
The correlations were taken into account by the second-order M0ller-Plesset perturbation 
theory (MP2) and by the coupled cluster method with single, double (and triple) cluster 
amplitudes, CCSD(T), for 19 external electrons of the E112H molecule. It is not clear 



from Ref. 
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why the correlations with the 6p shell of El 12 are considered but correlations 
with the 6s and 5/ shells are not whereas the small value of calculated dissociation energy, 
-D^ H =0.06eV, could be strongly influenced by the latters. The 6s and 6p shells are closely 
localized in space, whereas the 5/ and 6p shells have close orbital energies. The effect of the 
finite nuclear size was taken into account but the Breit effects and even the spin-orbit (SO) 
interaction were neglected (i.e., only scalar-relativistic calculations were made). It is clear 
that both effects are increased with the nuclear charge Z, therefore, these approximations 
can be inappropriate for SHE compounds even if they are justified for their lighter analogues. 

In Refs. y|]J, the RECP calculations of E112H+ (but not E112H) were carried out by 
the MP2 and CCSD(T) methods. In Ref. Q, complete active space SCF and multireference 
configuration interaction (MRCI) calculations of E112H + were also performed 4 ^. The SO, 
finite nuclear size and Breit effects were taken into account at the generation stage of the 
pseudopotential (PP) of Seth et alM. It should be noted that the parameters of this PP 
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were fitted with the help of the adjustment procedure based on the LS-coupling scheme 
(which was, in general, inappropriate for SHE's since the errors of these PP's were up to 
4 eV for valence energies in our test calculations accounting for the SO interaction non- 
perturbatively). Recently, a new PP for El 12 was generated by Seth et al. using the jj- 
coupling scheme^ 2 -, but we are not informed about molecular calculations with this PP. The 
SO interaction was taken into account in the calculations from Ref. 0. However, the author 
applied the 20-electron RECP of Nash et alM. The Breit effects were not considered at 
the generation stage of this RECP. It is not clear from Refs. which nuclear model was 
used there. In our test calculations 17 , the errors of this RECP in reproducing the results 
of the all-electron Dirac-Fock-Breit calculations with the Fermi nuclear model were up to 
1 eV for transition energies between low-lying states. Moreover, the basis set superposition 
errors (BSSE's) 23,24 were not estimated in the calculations from Ref. a. 

Methods and calculations. The GREC P 17 ' 18 ' 19 ' 25 , Fock-space relativistic CC-SD (RCC- 
SD) 26 ^ and spin-orbit direct configuration interaction (SODCI) 2 ^^ methods used for 
the present calculations are well described in literature. The gaussian expansions of our 
GRECP and (16, 21, 16, 12, 14)/[4, 6, 4, 2, 1] basis set for E112 are available at our website 
http:// www.qchem.pnpi.spb.ru/Basis/ In the SODCI calculation, the relativistic scheme 



of configuration selection was applied 30 . 

Two series of Fock-space RCC-SD calculations were performed for E112H with the 
GRECP. The ground state of the cation E112H + served as reference in the first series 
(denoted by RCC-SD- 1), and the Fock-space scheme was 

E112H+ ->• E112H, (1) 

with an electron added in the lowest unoccupied a orbital of E112H + . The second series 
(RCC-SD-2) started from the ground state of the anion E112H - as reference using the 
Fock-space scheme 

E112IT -»• E112H -> E112H+ (2) 

with electrons removed from the highest occupied a orbital of E112H~. Moreover, the RCC- 
SD calculations of El 12 (to calculate counterpoise corrections and D e ) were carried out where 
the 6d 10 7s 2 ground state of the El 12 atom was used as reference and the Fock-space scheme 
was 

E112->E112 + , (3) 



with an electron removed from the 6d or 7s shell. 

Our test atomic RCC calculations on E112 showed that at least 34 external electrons 
of the atom (occupying the 5/, 6s, 6p, 6d, 7s, . . . shells) should be correlated and the basis 
set should include up to z-type harmonics (/=6) in order to calculate the excitation and 
ionization energies with "chemical accuracy" . Nevertheless, we expect that the contributions 
of the core correlations will be less important for the molecule than for the atom as was in our 
similar calculation a 25 ' 31 ' 32 on Hg and HgH. This is, in particular, supported by a large orbital 
energy separation between hp and 5d shells in Hg (e[5p 3 / 2 ; 5g? 3 / 2 ] ~ —2.8; —0.65 a.u.) and by 
a comparable separation between 6p and 6d shells in E112 (e[6p3/2; 6^3/2] ~ —2.4; —0.56 a.u.) 
Therefore, only 13 external electrons for the E112H molecule (12 electrons for E112H + ) were 
correlated in the present calculations. The calculations with the larger number of correlated 
electrons when SO interaction is explicitly treated are rather time-consuming and suggested 
in future. 

In scalar-relativistic CC-SD calculations, we have also estimated (see Table HJ that 
correlations with the 6p shell of El 12 give relatively small contributions to the spectro- 
scopic constants in E112H and E112H + except for D e in E112H + . This cation disso- 
ciates to E112 + {Qd^ 2 6d^ 2 7sy 2 ) + H(ls 1 )^ in contrast to HgH + , which dissociates to 
Hg + (5o?3y 2 5c?5y 2 6s^ 2 ) + H(ls 1 ). The 6p shell is closely localized to the 6d shell, therefore, 
the correlations between these shells have to be important for transitions with an essential 
change in the occupation number for the 6c? shell. In principle, the D e value for E112H + can 
be easily corrected using our atomic RCC results for ionization potential of the 6g?5/2 sub- 
shell of E112. However, we observed large compensations between contributions accounting 
for correlations with core shells and for basis functions with high angular momenta. Thus, 
the above ionization potential from calculation with 12 correlated electrons in the basis in- 
cluding up to (yf-harmonics (/=4) differs from that with 52 electrons and I up to 8 only on 
+773 cm- 1 (0.10 eV). 

The calculations were carried out for 15 internuclear distances from 2.3 a.u. to 3.7 a.u. with 
interval of 0.1 a.u. The spectroscopic constants were calculated by the Dunham method in 
the Born-Oppenheimer approximation. All our RCC and SODCI results reported in Table HI 
were improved using counterpoise corrections (CPC's) 23 i 24 calculated for the El 12 6g? 10 7s 2 
state with a ghost H atom. CPC's calculated for the ground state of the H atom are about 
1 cm -1 , and are, therefore, ignored. 
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Results and discussion. Our results for the ground states of E112H and E112H + are 
collected in Table |U The corresponding result o 25 ' 32 for HgH and HgH + and the results 
of other group n 8 i 14 i 21 for E112H and E112H + are also presented for comparison. In the 
GRECP/RCC-SD-1 calculations one can observe the bond length contraction for E112H 
and E112H+ on 0.07 and 0.06 A with respect to HgH and HgH+. Detailed comparison of 
our results for HgH and HgH + with the results of other groups and the experimental data 



can be found in Ref. 



32. 



Our RCC-SD values for spectroscopic constants show considerable differences between 
two Fock-space schemes, GRECP/RCC-SD-1 and GRECP/RCC-SD-2. Such differences 
are caused by the truncation of the CC operator, they indicate significant contributions 
of the omitted higher-order (triple, etc.) cluster amplitudes (HOCA). HOCA influence on 
the total energies in each point of the potential curves were estimated with the help of 
the configuration interaction corrections on HOCA 34 calculated as differences in the total 
energies of the SODCI and RCC-SD-1 values. In these calculations, the same numbers of 
electrons were correlated as in the above RCC-SD case, but a reduced basis set, [4, 4, 3, 1] on 
E112 and [3,2] on H, was used because approaching the full configuration interaction limit 
in SODCI calculations becames too time-consuming for larger basis sets. 

Except for the dissociation limit, HOCA has small effect on E112H + since the cation 
is a closed shell system. It is well known that CC approach works particularly well for 
the closed-shell states which is confirmed by the comparison of the GRECP/12e- RCC-SD-1 
results with and without HOCA correction for E112H + . The change in D e is mainly due to 
different ionization potentials for the Qd^/2 electron of the E112 atom in the RCC-SD-1 and 
SODCI calculations. The GRECP/12e-RCC-SD-2 results show considerable dictinctions 
from the results corrected by HOCA because E112H + is calculated in the high Fock-space 
sector, (2,0), in which some lost of accuracy takes place. The HOCA contribution for 
E112H is important. Similar trend was observed in our GRECP/RCC calculations 32 on 
HgH when the effect of the triple cluster amplitudes was taken into account for 13 electrons 
that essentially improved the agreement with the experimental data. 

The differences between DK3/CCSD(T) 21 and our GRECP/RCC-SD-1 + HOCA results 
are small for E112H + , but are essentially larger for E112H. In particular, unlike our 
GRECP/13e- RCC-SD-1 + HOCA bond length, the DK3/19e-CCSD(T) value 21 for E112H is 
larger than the experimental dat a 1 ^ 3 ^ 3 ^ 3 ? for HgH. It worth to note that the DK3/CCSD(T) 
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calculations^ are scalar-relativistic whereas our RCC and SODCI calculations are performed 
with the spin-dependent GRECP. To check the effect of SO interaction, we have also car- 
ried out scalar-relativistic CC calculations with the spin-averaged GRECP part. The same 
basis, number of correlated electrons, Fock-space schemes, etc. were taken as in the RCC 
calculations. One can see from comparison of our CC-SD and RCC-SD results that the SO 
effect is small for E112H + (except for D e ) but is very essential for E112H. 

The R e and w e values for E112H+ by the PP/CCSD(T) method from Ref.Q (see the foot- 
note in TableHJ) differ from the DK3/18e-CCSD(T) results 21 by -0.017 A and +45 cm" 1 . The 
corresponding PP/CCSD(T)+SO values differ from our GRECP/12e-RCC-SD-l + HOCA 
results by -0.023 A and +126 cm -1 . The RECP/RCCSD(T) equilibrium distance for E112H+ 
from Ref. Q differs from our GRECP/RCC-SD-1 + HOCA result by +0.04 A and even more 
from the results of other group o 14 ' 21 . The difference between the RECP/RCCSD(T) equilib- 
rium distance calculated in Ref. for HgH + and the experimental datum is -0.04 A. Thus, 
one can observe the larger bond length for E112H + by 0.03 A in comparison with HgH + in 
the RECP/RCCSD(T) calculations 8 (in contrast to our results). 

One, however, can expect rather some increase in R e for E112H + and decrease for E112H 
when accounting for correlations with the innermore shells. In our scalar-relativistic calcu- 
lations, some small increase in D e was observed with enlarging the basis set. 

Conclusions. It is well known that properties of SHE's are somewhat different from 
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those of their lighter analogues due to very strong relativistic effects first of all (see Ref. 
and references). Therefore, even those approximations which work well for the lighter ana- 
logues (neglecting the SO interaction for £ states, the effective state of a considered atom in 
a molecule, the preferred valency, etc.) should not be used for SHE's without serious check- 
ing and analyzing. The calculated equilibrium distance, R e , in E112H is notably smaller 
than that in HgH. Therefore, one can also expect smaller bond lengths for the other El 12 
compounds in comparison with the Hg ones. There is a long-term discussion in scientific 
community whether E112 will behave like Hg or Rn. The ground state RnH and XeH 
molecules are not observed in the gas phase. Our calculations for the E112H molecule do 
not predict large dissociation energy, D e , but it is yet close to that of HgH. Therefore, we 
believe that the singly-valent ground state E112 compounds will rather resemble the Hg 
compounds than the noble gas patterns that is also supported by calculations of other El 12 
compounds 38,39,40 . 
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TABLE I: Spectroscopic constants of the ground states of the E112H molecule and the E112H + ion 
from two-component RCC-SD and scalar-relativistic CC-SD calculations with GRECP in the H 
(8, 4, 3)/ [4, 2,1] ANO and El 12 (16, 21, 16, 12, 14)/[4, 6, 4, 2, 1] basis set. Our corresponding results 
for HgH and HgH + and the results of other groups for E112H and E112H + are also presented for 



comparison. R e is in A, D e in eV, Y02 in 10 6 cm , other values in cm 



Molecula 


Method 


-Re We 


D e B e 




«e 


-Y 02 






Our calculations: 










HgH+ 


GRECP/12e-RCC-SD-l 


1.596 2037 


2.67 6.60 


39 


0.200 


279 


HgH+ 


GRECP/12e-RCC-SD(T)-l 


1.599 2013 


2.68 6.58 


41 


0.208 


282 


HgH+ 


Exp orimcnt a a " 1 ' 


1.594±0.000 2031±3 


(2.75±0.36) a 6.61±0.00 


44±3 


0.206±0.000 


285±0 



Our calculations: 



E112H+ 


GRECP/20e-CC-SD-l 


1.537 


2587 


4.60 


7.10 


46 


0.198 


215 


E112H+ 


GRECP/20e-CC-SD-2 


1.531 


2681 


4.46 


7.15 


35 


0.168 


205 


E112H+ 


GRECP/18e-CC-SD-l 


1.537 


2588 


4.61 


7.10 


17 


0.198 


215 


E112H+ 


GRECP/18e-CC-SD-2 


1.531 


2680 


4.46 


7.16 


35 


0.169 


205 


E112H+ 


GRECP/12e-CC-SD-l 


1.535 


2590 


4.96 


7.12 


47 


0.200 


216 


E112H+ 


GRECP/12e-CC-SD-2 


1.527 


2679 


4.75 


7.19 


37 


0.175 


208 


E112H+ 


GRECP/12e-RCC-SD-l 


1.537 


2569 


3.96 


7.11 


47 


0.201 


218 


E112H+ 


GRECP/12e-RCC-SD-2 


1.519 


2752 


3.80 


7.28 


15 


0.187 


204 


E112H+ GRECP/12e-RCC-SD-l - 


h HOCA 1.540 


2547 


4.35 


7.08 


45 


0.195 


220 






Other groups' 


calculations: 












E112H+ 


DK3/18e-CCSD21 


1.528 


2621 




7.18 








E112H+ 


DK3/18e-CCSD(T)2I 


1.532 


2595 




7.15 








E112H+ 


PP/CCSD(T)i£ b 


1.515 


2640 


5.15 




51 






E112H+ 


PP/MRCI+SOi4 b 


1.503 


2620 


3.86 










E112H+ 


PP/CCSD(T)+SOi£ b 


1.517 


2673 


4.09 




52 






E112H+ 


RECP/RCCSD(T)£ b 


1.583 




3.50 














Our calculations: 












HgH 


GRECP/13e-RCC-SD-l 


1.709 


1575 


0.35 


5.76 


56 


0.262 


312 


HgH 


GRECP/13e-RCC-SD(T)-l 


1.738 


1395 


0.41 


5.56 


83 


0.348 


363 


HgH 


Experiment 1D » 10 ' J0 


1.738±0.003 


1403±18 0.46±0.00 


5.57±0.02 


98±23 


0.337±0.067 


345±1 


HgH 


Experiment 37 


[1.766] c 


[1203] c 


0.46 


[5.39] c 






[395] c 






Our calculations: 












E112H 


GRECP/21e-CC-SD-l 


1.742 


1438 


-0.03 


5.53 


113 


0.409 


340 


E112H 


GRECP/21e-CC-SD-2 


1.801 


1104 


-0.02 


5.15 








E112H 


GRECP/19e-CC-SD-l 


1.741 


1439 


-0.03 


5.53 


113 


0.409 


340 


E112H 


GRECP/19e-CC-SD-2 


1.801 


1102 


-0.02 


5.15 








E112H 


GRECP/13e-CC-SD-l 


1.746 


1402 


-0.03 


5.50 


119 


0.429 


354 


E112H 


GRECP/13e-CC-SD-2 


1.808 


1038 


-0.05 


5.10 








E112H 


GRECP/13e-RCC-SD-l 


1.638 


1859 


0.36 


6.25 


95 


0.338 


288 


E112H 


GRECP/13e-RCC-SD-2 


1.663 


1649 


0.32 


6.06 


123 


0.425 


340 


E112H 


GRECP/13e-RCC-SD-l - 


h HOCA 1.662 


1800 


0.42 


6.07 


152 


0.385 


287 






Other groups' calculations: 












E112H 


DK3/19e-CCSD2i 


1.823 


991 


0.04 


5.05 








E112H 


DK3/19c-CCSD(T)2I 


1.829 


1007 


0.06 


5.02 









a Cited in Refs. 35 37 as uncertain. 

b Note that the RCCSD(T), RMRCI+SO, RCCSD(T)+SO values from Ref.[3and the CCSD(T) values from Ref.||are 
listed. The acronyms for these calculations (the last one is in the jj-coupling scheme, the other ones are scalar-relativistic 
where the second and third ones are corrected for the SO effects) are redefined in accordance with the other notations of the 
present paper. 

c Cited in Ref. l37l as corresponding to the zero vibrational level. 
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